TUM-HEP-347/99 
LMU 06/99 



Standard Model Confronting 



New Results for e'/e 

5—( ! 

S. Bosch 1 , A.J. Bums 1 , M. Gorbahn 1 , S. Jager 1 , 

< : 

M. Jamin 2 , M.E. Lautenbacher 1 and L. Silvestrini 1 

<n : 

1 Physik Department, Technische Universitat Miinchen, 
D-85748 Garching, Germany 

^ : 

O ■ Sektion Physik, Universitat Miinchen, 

G\ ■ 

Theresienstrasse 37, D-80333 Miinchen, Germany 

pu,: 

Qh! Abstract 

D . 

. We analyze the CP violating ratio e /e in the Standard Model in view of the new 

KTeV results. We review the present status of the most important non-perturbative 
parameters B^ 2 \ B^ 2 \ Bk and of the strange quark mass m s . We also briefly 
discuss the issues of final state interactions and renormalization scheme dependence. 
Updating the values of the CKM parameters, of m% and A^ and using Gaussian 
errors for the experimental input and flat distributions for the theoretical parameters 
we find e'/e substantially below the NA31 and KTeV data: e'/e = (7.7 ±|;§) • 1CT 4 
and e'/e = (5.2 j^'-jO'lO -4 m the NDR and HV renormalization schemes respectively. 
A simple scanning of all input parameters gives on the other hand 1.05 • 10 -4 < 
e'/e < 28.8 • 10~ 4 and 0.26 • 10~ 4 < e'/e < 22.0 • 10^ 4 respectively. Analyzing the 
dependence on various parameters we find that only for extreme values of B^ 2 \ 
Bf /2) and m s as well as suitable values of CKM parameters and A^-, the ratio 
e'/e can be made consistent with data. We analyze the impact of these data on the 
lower bounds for lmV t dV t *, Br(Ki — ► vr°z/P), Br{Ki — > 7r°e + e _ )di r and on tan/3 in 
the Two Higgs Doublet Model II. 



1 Introduction 



One of the most fascinating phenomena in particle physics is the violation of CP symmetry 
in weak interactions. In the Standard Model CP violation is supposed to originate in a 
single complex phase S in the charged current interactions of quarks This picture is 
consistent, within theoretical hadronic uncertainties, with CP violation in K° — K° mixing 
(indirect CP violation) discovered in K L — > 7T7t decays already in 1964 0] and described 
by the parameter e 0: 

e = (2.280 ±0.013) • l(T 3 exp(i$ E ), $ £ « j. (1.1) 

It is also consistent with the recent measurement of sin 2(3 from B — > ipKs at CDF H, 
although the large experimental error precludes any definite conclusion. 

It should be emphasized that the agreement of the Standard Model with the experi- 
mental value of e is non-trivial as | sin<5| < 1. Indeed in the Standard Model 

e = B K ImX t ■ F e (m t , ReA t ) exp(i7r/4) (1.2) 

where Bk is a non-pertubative parameter 0(1) and X t = V t dV t * with Vij being the elements 
of the CKM matrix |T], [|. The function F e results from well known box diagrams with 
W ± , t, c, u exchanges || and includes NLO QCD corrections 0,1. An explicit expression 
for F £ can be found in (10.42) of ||. It is an increasing function of the top quark mass 
m t and of ReA t . The QCD scale (Aj^g) dependence of F £ is very weak. 

Now, ImAt is an important quantity as it plays a central role in the phenomenology of 
CP violation in K decays and is furthermore closely related to the Jarlskog parameter Jqp 



10|1, the invariant measure of CP violation in the Standard Model: Jqp = Ayl — A 2 ImA 4 



with A = 0.221 denoting one of the Wolfenstein parameters [0. To an excellent approx- 
imation one has 

ImA t = iKbHK&l sin5. (1.3) 



As can be inferred from (|1.2j ) and (|1.3| ) only for sufficiently large values of \V u b\ } |V^,|, 
m t and Bk can e in ( |1.2|) and consequently the indirect CP violation in the Standard 
Model be consistent with the one observed experimentally. It turns out that using the 
known values of \V u b\, \V c b\, m t (see Section 3) and taking Bk = 0.80 ±0.15 in accordance 
with lattice and large-N calculations (see Section 2), the experimental value of e can be 
reproduced in the Standard Model provided sin 5 > 0.69. This determination of sin S 
includes constraints from B® s — B^ s mixings. We also find 

1.04- 10" 4 < ImA t < 1-63- 10~ 4 . (1.4) 
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It should be noticed that sin 8 = 0(1) and that the extracted range for lm\ t is not 
far from the upper limit of 1.73 • 10~ 4 following from the unitarity of the CKM matrix. 
It should also be emphasized that the large top quark mass plays an important role in 
obtaining the experimental value for e. Had m t been substantially lower than it is, the 
theoretical value of e would be below the experimental one. 

While indirect CP violation in Kl — > tttt reflects the fact that the mass eigenstates in 
the K° — K° system are not CP eigenstates, the so-called direct CP violation is realized 
via direct transitions between states of different CP parities: CP violation in the decay 
amplitude. In Kl —* vr7r decays this type of CP violation is characterized by the parameter 
e'. In the Standard Model one has 



£ 



-=lm\ t -F e/ (m t ,A^,m s ,Bi l/2 \Bi 3/2 \n ri+ril ) (1.5) 
where the function F £ i results from the calculation of QCD penguin and electroweak 

(1/2) (3/2) 

penguin diagrams. Here B\' 1 and Bg 1 ' are non-perturbative parameters related to the 
dominant QCD penguin and electroweak penguin contributions respectively, is the 
QCD scale and VL V+V i represents isospin breaking effects. 

The expression ( |1.5| ) has been obtained by calculating e' and dividing it by the ex- 
perimental value of e in ( |1 . 1| ) in order to be able to compare with the experimental value 
of e'/e. This procedure exhibits the nature of e' which representing direct CP violation 
is proportional to lm\ t . However, one could also proceed differently and ignoring the 
constraint ( |L1|) calculate e'/e fully in theory. In this case ( |1.5|) is replaced by 

^ _ ^(m t ,A|, mg ,gj 1/2 Ug ,/2) ,n^) 



e B K F £ (m t ,ReX t ) 

where F £ / = |e exp |.F e / is independent of e. One should notice that ImAt cancelled out in 
£'/£ calculated in this manner and e'/e is actually a function of KeX t and not of lm\ t . 
However, once the constraint ( |1.1| ) has been taken into account ( |1.6| ) reduces to (|L5p. We 
will return to this point in Section 3. 

There is a long history of calculations of e'/e in the Standard Model. The first cal- 
culation of e'/e for m t <C without the inclusion of renormalization group effects can 
be found in [ l2"| . Renormalization group effects in the leading logarithmic approximation 
have been first presented in fl3| . For m t <^ M\y only QCD penguins play a substantial 



role. First extensive phenomenological analyses in this approximation can be found in [14 



Over the eighties these calculations were refined through the inclusion of QED penguin 
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effects for m t <C Mw |T3|, 0, |T7[, the inclusion of isospin breaking in the quark masses 
fl6| , |17], [181 ) and through improved estimates of hadronic matrix elements in the frame- 



work of the 1/N approach ||19|| . This era of e'/e culminated in the analyses in p(| |21 



where QCD penguins, electroweak penguins (7 and Z° penguins) and the relevant box 
diagrams were included for arbitrary top quark masses. The strong cancellation between 
QCD penguins and electroweak penguins for m t > 150 GeV found in these papers was 



confirmed by other authors [22]. 

During the nineties considerable progress has been made by calculating complete NLO 
corrections to e' ||23||-||27||. Together with the NLO corrections to e and B°—B° mixing ]7|,[3, 
28 1 , this allowed a complete NLO analysis of e'/e including constraints from the observed 
indirect CP violation (e) and B® s — Bj s mixings (AM^ S ). The improved determination of 
the Vub and V c b elements of the CKM matrix, the improved estimates of hadronic matrix 
elements using the lattice approach as well as other non-perturbative approaches and in 
particular the determination of the top quark mass m t had of course also an important 
impact on e'/e. 

In a crude approximation (not to be used for any serious analysis) 



F F , « 13 ■ 



HOMeV 2 
m s (2 GeV) 



d(V 2 ) (-1 O \ OA p(3/2) ( m t 



2.5 



MS 



340 MeV 



(1.7) 

where VL v+rj i w 0.25. This formula exhibits very clearly the dominant uncertainties in F e i 
which reside in the values of m s , B^ 2 \ B^ 2 \ and £l v +rf- Because of the accurate 
value m t (m t ) = 165 ± 5 GeV, the uncertainty in e'/e due to the top quark mass amounts 
only to a few percent. A more accurate formula for F e i will be given in Section 2. 

A comparison of the formulae ( [L.2|) and ( p..5|) reveals that the analysis of e is theoreti- 
cally cleaner. Indeed, e depends on a single non-perturbative parameter Br, whereas e'/e 
is a sensitive function of B^ 2 \ B^ 2 \ m B , A^ and VL V+V i . Moreover, the partial can- 
cellation between QCD penguin (Bq ) and electroweak penguin {B% ) contributions 
requires accurate values of B^ 2 ^ and B^ 2 ^ for an acceptable estimate of e'/e. 

Until recently the experimental situation on e'/ e was rather unclear. While the result 
of the NA31 collaboration at CERN with Re(e'/e) = (23.0 ± 6.5) ■ 10~ 4 |2| clearly 
indicated direct CP violation, the value of E731 at Fermilab, Re(e'/e) = (7.4 ± 5.9) • 10~ 4 
SOD , was compatible with superweak theories |ST| in which e'/e = 0. This controversy is 
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now settled with the very recent measurement by KTeV at Fermilab [52| 

Re(-) = (28.0 ±4.1) ■ 10~ 4 (1.8) 

which together with the NA31 result confidently establishes direct CP violation in nature. 
The grand average including NA31, E731 and KTeV results reads 

pi 

Re(-) = (21.8 ±3.0) • 10" 4 (1.9) 

very close to the NA31 result but with a smaller error. The error should be further 
reduced once the first data from NA48 collaboration at CERN are available and complete 
data from both collaborations have been analyzed. It is also of great interest to see what 
value for e'/e will be measured by KLOE at Frascati, which uses a different experimental 
technique than KTeV and NA48. 

Does the direct CP violation observed in — ► titt decays agree with the Standard 
Model expectations? Before entering the details let us take a set of "central" values for the 
parameters entering F £ r. Together with B K = 0.80, m t {m t ) = 165 GeV, \V u b\ = 3.56- 10~ 3 
and \Vcbl — 0.040 needed for the e- analysis we set 

^ 1/2) = 1.0, ^ 3/2) =0.8, m s (2 GeV) = 110 MeV, Q v+rl , = 0.25 (1.10) 

and AJJL = 340 MeV. Using the formula (|2~38|) for F £ /, we find F e , = 5.2. On the other 
hand the e-analysis gives ImAt = 1.34 ■ 10~ 4 . Consequently 

central 



E 



7.0 -10" 4 (1.11) 



well below the experimental findings in ( |1.9| ). 

Equivalently, with F e > = 5.2, the experimental value in (|1.9|) implies ImA( = (4.2 ± 
0.6) ■ 10" 4 which lies outside the range (|1.4|) extracted from the standard analysis of the 
unitarity triangle. Moreover it violates the upper bound ImA t = 1.73- 10 -4 following from 
the unitarity of the CKM matrix. 

The fact that for central values of the input parameters the size of e 1 / e in the Standard 
Model is well below the NA31 value of (23.0±6.5)-10 -4 has been known for some time. The 

(1/2) (3/2*) 

extensive NLO analyses with lattice and large-N estimates of Bq ' w 1 and B\ ' ' « 1 
performed first in |25], ^] and after the top discovery in [Q-[^5] have found e'/e in 
the ball park of (3 - 7) ■ 10~ 4 for m s (2 GeV) w 130 MeV. On the other hand it has 
been stressed repeatedly in || |34| that for extreme values of B^ 2 \ B^ 2 ^ and m s still 
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consistent with lattice, QCD sum rules and large-N estimates as well as sufficiently high 
values of lm\ t and A^, a ratio e'/e as high as (2 — 3) • 10~ 3 could be obtained within the 
Standard Model. Yet, it has also been admitted that such simultaneously extreme values 
of all input parameters and consequently values of e'/e close to the NA31 result are rather 
improbable in the Standard Model. Different conclusions have been reached in , where 
values (1 — 2) ■ 1CT 3 for e'/e can be found. Also the Trieste group |38|, which calculated 
the parameters B^ 2 ^ and B^ 2 ^ in the chiral quark model, found e'/e = (1.7 ± 1.4) ■ 10~ 3 . 
On the other hand using an effective chiral lagrangian approach, the authors in |39[ found 
e'/e consistent with zero. 



The purpose of the present paper is to update the analyses in || |34| and to confront 



the Standard Model estimates of e'/e with the experimental findings in (|1.9| ). Other very 



recent discussions of e'/e can be found in |40||-||42||. We will comment on them below. In 



the present paper we address in particular the following questions: 



What is the maximal value of e'/ e in the Standard Model consistent with the usual 

MS 



analysis of the unitarity triangle as a function of B^ 2 \ B^ 2 \ m a and A^- ? 



What is the lowest value of B^ 2 ^ as a function of B^ 2 ^ for fixed values of m s and 
for which the Standard Model is simultaneously compatible with ( |1.9| ) and the 
analysis of the unitarity triangle? 

What is the sensitivity of the analysis of e'/e to the values of £l, n+ , n i and Bj<p. 

What is the impact of the experimental value for e'/e on ImA(, on the usual analysis 
of the unitarity triangle and in particular on Standard Model expectations for the 
rare decays Kl — > -n^vv and Kl — > 7r°e + e~ in which direct CP violation plays an 
important role? 



What are the general implications of ( |1.9| ) for physics beyond the Standard Model? 



In particular, what is the impact on the allowed range in the space (MH,tan/3) in 
the so called two Higgs doublet model II (2HDMII) 



While addressing these questions we would like to emphasize that it is by no means 
the purpose of our paper to fit Bq 2 ', B^ 2 \ m s , A^, VL V+V i and Bk in order to make the 
Standard Model compatible simultaneously with experimental values on e'/e, e and the 
analysis of the unitarity triangle. Such an approach would be against the whole philosophy 
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of searching for new physics with the help of loop induced transitions as represented by 
e'/e and e. Moreover it should be kept in mind that: 

• B^ 2 \ B^ 2 ^ and B K , in spite of carrying the names of non-perturbative parameters, 
are really not parameters of the Standard Model as they can be calculated by means 
of non-perturbative methods in QCD. The same applies to Q v+r) >. 

• m s , A^||j, m t , \Vcb\ and \V u t,\ are parameters of the Standard Model but there are 
better places than e'/e to determine them. In particular the usual determinations 
of these parameters can only marginally be affected by physics beyond the Standard 
Model, which is not necessarily the case for e and e'/e. 

Consequently, the only parameter to be fitted by direct CP violation is sin 5 or ImAt. 
The numerical analysis of e'/e as a function of B^ 2 \ B^ 2 \ m s , A^, Q v+T f and Bp- 
should only give a global picture for which ranges of parameters the presence of new 
physics in e'/ e and e should be expected. 

Our paper is organized as follows. In Section 2 we recall briefly the basic formulae 
for e'/e in the Standard Model. We also review the existing methods for estimating 
hadronic matrix elements of relevant local operators and we present a rather accurate 
analytic formula for F £ >. In Section 3 we address several of the questions listed above. In 
Section 4 we discuss briefly general implications for physics beyond the Standard Model. 
In particular we investigate the lower bound on tan (3 as a function of the charged Higgs 
mass in the 2HDMII. Conclusions and outlook are given in Section 5. 

2 Basic Formulae 
2.1 Formulae for e'/e 

The parameter e' is given in terms of the isospin amplitudes Ai as follows 

e' = -±=lm exp(^), ®e> = \ + h - 5 , (2.1) 

where 5j are final state interaction phases. Then, the basic formula for e'/e is given by 

^=ImVF £ ,, (2.2) 

where 

F e ,= [P(V2) _ P (3/2)] exp( (2.3) 
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with 



pd/2) = rX)yi<Oi>o(l-JW), (2-4) 
p(m = Lj2yi(Qi} 2 . (2.5) 



Here 
Since 

$ = $ e , - $ e « 0, (2.7) 

F £ / and e'/e are real to an excellent approximation. 
The operators Qi are given explicitly as follows: 
Current Current : 

Ql = {s a Up)v-A (U/3d a )v-A Q2 = (SU) V - A (ud)y_ A (2.8) 

Q CD— Penguins : 

Q3 = {sd) V -A {Q0)v-a Qa = {s a dp)v- A Y {qpq a )v-A (2.9) 

q=u,d,s q=u,d,s 

Q5 = {sd)v- A Y (^q)v+a Q% = {s a dp)v- A Y iWa)v+A (2.10) 

q=u,d,s q=u,d,s 

Elect roweak— Penguins : 

3 3 

Q7 = 2 Mv-A J2 e q(Ql)v+A Qs = ^(S a df3)v-A Y e q(W<*)v+A (2.11) 
q=u,d,s q=u,d,s 

3 3 
Q9 = 2 (sd)v-A Y e q(q<l)v-A Q10 = 2 (s a dp) V -A Zq(Wo)v-A- (2.12) 

" q=u,d,s " q=u,d,s 

Here, a, /3 are colour indices and e q denotes the electric quark charges reflecting the 
electroweak origin of Q 7 , . . . , Q w . 

The Wilson coefficient functions yi(/J,) were calculated including the complete next- 
to- leading order (NLO) corrections in p3|-|p7|. The details of these calculations can be 



found there and in the review fP| . Their numerical values for corresponding to 



a|^(Mz) = 0.119 ± 0.003 and two renormalization schemes (NDR and HV) are given 
in table [IJ. There we also give the coefficients z\$ relevant for the discussion of hadronic 
matrix elements. 
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Table 1: AS" = 1 Wilson coefficients at /i = m c = 1.3 GeV for m t = 165 GeV and / = 3 
effective flavours. y 1 = y 2 = 0. 





A^ = 290MeV 


A^ = 340MeV 


Ag^ = 390 MeV 


Scheme 


NDR 


HV 


NDR 


HV 


NDR 


HV 


Z\ 


-0.393 


-0.477 


-0.425 


-0.521 


-0.458 


-0.570 


Z2 


1.201 


1.256 


1.222 


1.286 


1.244 


1.320 


2/3 


0.027 


0.030 


0.030 


0.034 


0.033 


0.038 


2/4 


-0.054 


-0.056 


-0.059 


-0.061 


-0.064 


-0.067 


2/5 


0.006 


0.015 


0.005 


0.016 


0.003 


0.017 


2/6 


-0.082 


-0.074 


-0.092 


-0.083 


-0.105 


-0.093 


Vi/at 


-0.038 


-0.037 


-0.037 


-0.036 


-0.037 


-0.034 




0.118 


0.127 


0.134 


0.143 


0.152 


0.161 




-1.410 


-1.410 


-1.437 


-1.437 


-1.466 


-1.466 


W« 


0.496 


0.502 


0.539 


0.546 


0.585 


0.593 



It is customary in phenomenological applications to take ReAo and u from experiment, 



i.e. 



ReAn = 3.33 • 10~ 7 GeV, 



to 



0.045, 



;2.13) 



where the last relation reflects the so-called AI =1/2 rule. This strategy avoids to a large 
extent the hadronic uncertainties in the real parts of the isospin amplitudes Aj. In order 
to be consistent the constraint ( |2.13| ) should also be incorporated in the matrix elements 
(Qi)i necessary for the evaluation of e'/e. This in fact has has been done in 123] and we 
will return to this approach briefly below. Studies of the A J = 1/2 rule can be found in 



45, 46, 47 



The sum in ( |2.4j ) and ( j2.5| ) runs over all contributing operators. p( 3 / 2 ) is fully dom- 
inated by electroweak penguin contributions. P^ 1 / 2 ) on the other hand is governed by 
QCD penguin contributions which are suppressed by isospin breaking in the quark masses 
{m u ^ m d ). The latter effect is described by 



1 (ImA 



2J1.B. 



r)+rf 



10 



For Vt v+rj i we will first set 



n 



v+v' 



ImA 



0.25 



(2.14) 
(2.15) 
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which is in the ball park of the values obtained in the 1/N approach |T7| and in chiral 
perturbation theory |0| pl|. Vt rq+V i is independent of m%. We will investigate the sensitivity 
of e'/e to fln+rj' in Section 3. 



2.2 Hadronic Matrix Elements 

The main source of uncertainty in the calculation of e' / e are the hadronic matrix elements 
(Qi)i- They generally depend on the renormalization scale \x and on the scheme used to 
renormalize the operators Qi. These two dependences are canceled by those present in the 
Wilson coefficients yi(fi) so that the resulting physical e'/e does not (in principle) depend 
on \i and on the renormalization scheme of the operators. Unfortunately, the accuracy of 
the present non-perturbative methods used to evalutate (Q%)i is not sufficient to have the 
H and scheme dependences of (Qi)i fully under control. We believe that this situation 
will change once the lattice calculations and QCD sum rule calculations improve. A brief 
review of the existing methods including most recent developments will be given below. 

In view of this situation it has been suggested in to determine as many matrix 
elements {Qi)i as possible from the leading CP conserving K — > ttti decays, for which 
the experimental data is summarized in (|2.13|) . To this end it turned out to be very 
convenient to determine (Qi)i in the three-flavour effective theory at a scale \x ~ m c . 
With this choice of fi the operators Q\ 2 i being present only for \i > m c , are integrated 
out and the contribution of penguin operators to ReAi turns out to be very small. Un- 
fortunately, since the charm mass is not much larger than the scale Mr of the process 
we are studying, the matching procedure between the four- and three-flavour effective 
theories contains an ambiguity related to the choice of external momenta in the matching 
2"5|. Furthermore, as pointed out in P3[, there is an ambiguity due to the contribu- 



tion of higher dimensional operators which are unsuppressed for /x pa m c . However, all 
these ambiguities are of 0(a s ) and one can easily verify that their possible contribution 
to ReAi is at the level of a few percent at most. Consequently, they have only a minor 
impact on our determination of (Qi)i at jj, — m c from ReAi. Using the renormalization 
group evolution one can then find (Qi)i at any other scale fi ^ m c . The details of this 
procedure can be found in |2f| . 

As we will see below this method allows to determine only the matrix elements of 
the (V — A) <g> (V — A) operators. For the central value of lm\ t these operators give a 
negative contribution to e'/e of about —2.5 • 10~ 4 . This shows that these operators are 
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only relevant if e'/e is below 1 • 10 -3 . Unfortunately the matrix elements of the dominant 
(V — A) <g> (V + A) operators cannot be determined by the CP conserving data and one 
has to use non-perturbative methods to estimate them. 

Before giving the results for (Qi)i in our approach we would like to emphasize why 
it is reasonable to extract hadronic parameters from Re Aj, while this would not be the 
case for 1mA j, which govern e'/e. The point is that KeAj, in contrast to ImAj, are not 
expected to be affected by new physics contributions. 

It is customary to express the matrix elements (Qi)i in terms of non-perturbative 
parameters B^ and B^ 2 ^ as follows: 



(vac) 



{Qi) 2 = B? /2) (Qi} { 2 



(vac) 



(2.16) 



The label "vac" stands for the vacuum insertion estimate of the hadronic matrix elements 



in question for which B\ 



(1/2) 



B 



(3/2) 



Then the approach in gives at /i = m c : 

0.187GeV 3 



(Qi(m c )) = 

(<?20c))o -- 

{Qz{m c )) -- 
((5 4 (m c )) = 
(Q 5 (m c )) = 

(Qe(m c )) ~- 

(Q 7 (m c )) ~- 

{Q8(m c )) ~- 

(g 9 (m c )) = 

(Q 10 (m c )) ~- 

(Qi{m c )) 2 

(Qih 
(Q 7 (m c )) 2 



z 2 [m c 



(Q2{m c )) 



*XBi 1/2) 
9 2 



(rn c ) 



-XB. 



(1/2), 



m. 



(Q3(m c )) + (Q 2 (m c )) - (Qi{m c )) , 
^J 1/2) (m c )(Q^K)) , 

-4W 



m s {m c ) + m d (m c ) 
X 



(1/2), 



K 



(Q 6 (m c ))o(/« + i; 



2 



4 1/2) K), 



4 1/2) K), 



-(Q 6 (m c )) (fi; + 1) - — 

3 1 

-(Qi(m c )) - -(Q 3 (m c ))o, 

(Q2(m c )) + -(Qi(m c )) - ^(Q 3 (wO)o 



8.44 • 1(T 3 GeV 3 



(<5 2 (^c))2 , . 

z + (m c ) 
0, i = 3,..., 6, 



6^<QaK)>o + ;l 



(3/2), 



m c ) 



2.17) 

2.18) 

2.19) 
2.20) 
2.21) 

2.22) 

2.23) 
2.24) 
2.25) 
2.26) 

2.27) 
2.28) 
2.29) 
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(QsOc))2 



K 



(Qe(m c )) + 



2V5^ WU 6 



4 3/2) (m c ), (2-30) 



where 



(Qg(m c )) 2 = (Qio(m c )) 2 =-(Qi(m c )> 2j (2.31) 



K X = J 3 -F 7r (ml-ml) J (2.32) 



and 



Fk~ F-k V 2 

/7T7 \\ (Q6( m c))o , n qq\ 

(Q6(m c )) = (1/2) — - , ^+ = *i + Z2- (2.33) 
-Be ( m c) 

The equality of the matrix elements in ( |2.27| ) follows from isospin symmetry of strong 
interactions. Finally, by making the very plausible assumption, valid in known non- 
perturbative approaches, that (Q-(m c ))o > (Q + (m c ))o > 0, where Q± = (Q2 ± Qi)/2, 
B 2 (m c ) can be determined as well. This gives for = 340 MeV 

Bi% 2) DR (m c ) = 6.5 ± 1.0, Bi% 2 l{m c ) = 6.1 ± 1.0 . (2.34) 

The actual numerical values used for uik, Tn n , F K , F n are collected in the appendix of 



4||. In particular = 131 MeV. 

It should be noted that this method allows to determine not only the size but also 
the renormalization scheme dependence of those matrix elements which can be fixed in 
this manner. This dependence enters through z 12 (m c ) and the scheme dependence of 
B^ 2 ^ (m c ). In obtaining the results above one also uses operator relations valid for /x < m c 
which allow to express Q 4 , Q 9 and Q w in terms of Qi, Q 2 and Q 3 . Theoretical issues 
related to these relations in the presence of NLO QCD corrections and the case of matrix 
elements for \i > m c are discussed in detail in |25[ . 



In order to proceed further one has to specify the remaining Bi parameters in the 
formulae above. As the numerical analysis in ]23| shows e'/e is only weakly sensitive 



to the values of the parameters B^ 2 \ B^ /2 \ B^ /2 \ B { s 1/2) and B? /2) as long as their 



absolute values are not substantially larger than 1. As in [25] our strategy is to set 



bS!(?K) = 1, Bi 1/2 \m c ) = Bi 1/2 \m c ), B? /2) (m c ) = B$ l2 \m c ) (2.35) 
and to treat Bq (m c ) and B^ 2 \m c ) as free parameters. 



The approach in |25| allows then in a good approximation to express e'/e or equiv- 
alently F e > in terms of A^, m t , m s and the two non-perturbative parameters B^ 2 ^ = 
Bq (m c ) and B^ 2 ^ = B% (m c ) which cannot be fixed by the CP conserving data. 
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2.3 The Issue of Final State Interactions 



In ( |2.1|) and (|2.7 ) the strong phases 5$ ~ 37° and 52 ~ —7° are taken from experiment. 



They can also be calculated from NLO chiral perturbation for titt scattering [37]. How- 



ever, generally non-perturbative approaches to hadronic matrix elements are unable to 



reproduce them at present. As 5i are factored out in (|2.1|) , in non-perturbative calcu- 
lations in which some final state interactions are present in (Qi)i one should make the 
following replacements in ( [2.4| ) and ( |2.5| ): 

{Qi)i - (2.36) 
(cosd/jth 

in order to avoid double counting of final state interaction phases. Here (cos<5/) t h is 
obtained in a given non-perturbative calculation. In leading large-N calculations and 
in quenched lattice calculations the phases 5j vanish and this replacement is ineffective. 
When loop corrections in the large-N approach [|7j, 4~3] and in the chiral quark model 



38fl are included an absorptive part and related non- vanishing phases are generated. Yet, 



in most calculations the phases are substantially smaller than found in experiment. For 
instance in the chiral quark model (cos<5o)th ~ 0.94 to be compared to the experimental 
value (cos<5o)exp ~ 0.8. Even smaller phases are found in |L9], gTJ, 

The above point has been first discussed by the Trieste group who suggested that 



in models in which at least the real part of {Qi)i can be calculated reliably, one should 
make the following replacements in (|2.4| ) and 



(Qi)i - (2.37) 

(cOSf)/) oxp 

where this time the experimental value of Si enters the denominator. As (cos5o) CX p ~ 0.8 
and (cosi^exp ~ 1 this modification enhances P^ 1 / 2 ) by 25% leaving p( 3 / 2 ) unchanged. 
The same procedure has been adopted in ^7j. To our knowledge there is no method for 
hadronic matrix elements which can provide 5q ~ 37° and consequently the replacement 
( p.37| ) may lead to an overestimate of the matrix elements. 



As in our paper the matrix elements of (V — A) ® (V — A) operators are extracted from 
the data, the replacements in (|2.36|) and ( j2.37| ) are ineffective for the determination of 
the corresponding contributions. They merely change the definition of the Bi parameters 
in the matrix elements of (V — A) <g> (V — A) operators. The situation is different with 
the matrix elements of (V — A) <g> (V + A) operators which are taken from theory. Yet in 
view of the remarks made above, in our analysis we will use exclusively (|2.4j) and (|2.5|) 
including possible effects of this sort in the uncertainties in Bq'^ and B^ 2 \ 
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2.4 An Analytic Formula for e'/e 



As shown in [49j], it is possible to cast the formal expressions for e'/e in (|2.2j )— (|2^ ) into an 
analytic formula which exhibits the mt dependence together with the dependence on m s , 
A^, B^ 2 ^ and B^ 2 ^ . To this end the approach for hadronic matrix elements presented 
above is used and Q v +-q' is set to 0.25. The analytic formula given below, while being 
rather accurate, exhibits various features which are not transparent in a pure numerical 
analysis. It can be used in phenomenological applications if one is satisfied with a few 
percent accuracy. Needless to say, in our numerical analysis in Section 3 we have used 
exact expressions. 

In this formulation the function F e > is given simply as follows (xt = mf/M^)' 

F E , =P + Px X {x t ) + P Y Y (x t ) + P z Z (x t ) + P E E (x t ). (2.38) 
Exact expressions for the mt-dependent functions in ( |2.38| ) can be found for instance in 



H |44|. In the range 150 GeV < m t < 180 GeV one has to an accuracy much better than 
1% 

m t 



Xo(x t ) = 1.51 



Z (x t ) = 0.66 



1.13 



165 GeV, 

m t \ L90 
165 GeV J ' 



Y (xA 



E (x t ) 



-- 0.96 
0.27 



mt 



1.55 



In our numerical analysis we use exact expressions. 

The coefficients Pj are given in terms of B^ 2 ^ = B^ 2 ^ 
mJm r ) as follows: 



165 GeV 

mt 
165 GeV, 

m c ), B^ 



-1.08 



B 



(3/2), 



JO) , (6) p (8) R 



where 



p - d(1/2) 



137 MeV 



1 2 



m s (m c ) + m d (m c ) 



i?s = B, 



(3/2) 



137 MeV 



m s (m c ) + m d (m c ) 



(2.39) 
(2.40) 

m c ) and 
(2.41) 

(2.42) 



The Pi are renormalization scale and scheme independent. They depend, however, on 
A^. In table |2] we give the numerical values of rf \ rf^ and r,- for different values of 
1 1 '■ - m r in the NDR renormalization scheme. This table differs from the ones 



presented in || [34]] in the values of A^ and the central value of m s (m c ) in R s which 



MS 



has been lowered from 150 MeV to 130 MeV. The coefficients rf\ rf^ and rf^ depend 
only very weakly on m s (m c ) as the dominant m s dependence has been factored out. The 
numbers given in table ^| correspond exactly to m s (m c ) = 130 MeV. However, even for 
m s (m c ) ~ 100 MeV or m s (m c ) ~ 160 MeV, the analytic expressions given here reproduce 
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the numerical calculations of e'/e given in Section 3 to better than 4%. For different scales 
\i the numerical values in the tables change without modifying the values of the Pj's as 



it should be. The values of B^ 2 ^ and B^ 1 *' should also be modified, in principle, but 



,(3/2) 



as a detailed numerical analysis in 2o| showed, it is a good approximation to keep them 
/^-independent for 1 GeV < fi < 2 GeV. We will return to this point below. 



Table 2: Coefficients in the formula ( |2.41| ) for various in the NDR scheme. The last 
row gives the r coefficients in the HV scheme. 



A^ = 290MeV 



A^ = 340MeV 



A||= 390 MeV 



(0) 



,(6) 



(0) 



,(6) 



.(8) 



(0) 



(6) 



,(8) 





X 
Y 
Z 

Eq 





-2.771 
0.532 
0.396 
0.354 
0.182 
-2.749 



9.779 
0.017 
0.072 
-0.013 
-1.144 
8.596 



1.429 



-9.404 
0.411 
1.050 



-2.811 
0.518 
0.381 
0.409 
0.167 

-2.788 



11.127 
0.021 
0.079 
-0.015 
-1.254 
9.638 



1.267 



-10.230 
0.461 
0.871 



-2.849 
0.506 
0.367 
0.470 
0.153 

-2.825 



12.691 
0.024 
0.087 
-0.017 
-1.375 
10.813 



1.081 



-11.164 
0.517 
0.669 



The inspection of table ^ shows that the terms involving r$ and r z dominate the 
ratio e'/e. Moreover, the function Z (x t ) representing a gauge invariant combination of 
Z°- and 7-penguins grows rapidly with m t and due to < these contributions suppress 
e'/e strongly for large m t |20|, |2T|j . 



2.5 Renormalization Scheme Dependence 

Concerning the renormalization scheme dependence only the coefficients j-q , Tq and 
are scheme dependent at the NLO level. Their values in the HV scheme are given in 
the last row of table [| We note that the parameter Tq is essentially the same in both 
schemes as the dominant scheme independent contributions to have been determined 
by the data on ReAj. Since Pq must be scheme independent and and are scheme 

(1/2) (3/2^ 

dependent, we conclude that Bq and B$ must be scheme dependent. Indeed the 
matrix elements in the NDR and HV schemes are related by a finite renormalization which 



can be found in equation (3.7) of ||25|| . Using this equation together with the approach 
to matrix elements presented above, we find approximate relations between the values of 
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(fi 6 , fig ) in the NDR scheme and the corresponding values in the HV scheme: 

(^ 1/2) )hv w 1.2(B? /2) )ndr, (Sf /2) )hv « 1.2(sf /2) )ndr- (2.43) 

One can check that the scheme dependence of (B^ 2 \ fig 3 ^) cancels to a very good 
approximation the one of and so that Pq is scheme independent. 

On the other hand the coefficients r*, i — X,Y, Z, E are scheme independent at NLO. 
This is related to the fact that the m t dependence in e'/e enters first at the NLO level 
and consequently all coefficients r« in front of the mt dependent functions must be scheme 

(1/2) (3/2) 

independent. Strictly speaking then the scheme dependence of fig and fig inserted 
into Pi with % ^ is really a part of higher order contributions to e'/e and should 
be dropped at the NLO level. Formally this can be done by not performing the finite 
renormalization when going from the NDR to the HV scheme. Then the coefficients Pj 
with % ^ are clearly scheme independent. 

In practice the situation is more complicated. The present non-perturbative methods 

(1/2) (3/2) 

used to evaluate fi 6 and fig like the large-N approach are not sensitive to the 
renormalization scheme dependence and we do not know which renormalization scheme 
the resulting values for these parameters correspond to. Lattice calculations, QCD sum 
rule calculations and the chiral quark model can in principle give us the scheme dependence 

(1/2) (3/2) 

of fi^ ; and B^' J but the accuracy of these methods must improve before they could 
be useful in this respect. 

In view of this situation our strategy will be to use the same values for fig 1//2 ' ) and 

(3 /2) 

fig in the NDR and HV schemes. This will introduce a scheme dependence in Pq and 
consequently in e'/ e but will teach us something about the uncertainty in e' / e due to the 
poor sensitivity of present methods to renormalization scheme dependence. 

It should also be noted that even if we knew the scheme dependence of fi 6 1//2 ' ) and 

(3/2) 

fi 8 without the ability of separating a scheme independent part in these parameters, 
the resulting e'/e would be scheme dependent at the NLO level. This time the scheme 
dependence would enter through the scheme dependence of Pi with % ^ 0. The latter 
scheme dependence could only be reduced by including the next order of perturbation 
theory in the Wilson coefficients: a formidable task. We should also stress [25[ that the 
scheme dependences discussed here apply not only to QCD corrections but also to QED 
corrections. That is QED corrections to the matrix elements of operators have to be also 
known. 

For similar reasons the NLO analysis of e' / e is still insensitive to the precise definition 
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of rat. In view of the fact that the NLO calculations needed to extract ImAt have been 
performed with ra t = ra t (ra t ) we will also use this definition in calculating F e i. 



2.6 Status of the Strange Quark Mass 

At this point it seems appropriate to summarize the present status of the value of the 
strange quark mass. Since different methods provide ra s at different values of \i we give 
in table |3| a dictionary between the ra s values at \i = 1 GeV, /i = ra c = 1.3 GeV and 
H = 2 GeV. 

In the case of quenched lattice QCD the present status has been summarized recently 
by Kenway [|5(J. Averaging the results presented by him at LATTICE 98, we obtain 



ra s (2 GeV) = (120 ± 20) MeV. It is expected that unquenching will lower this value but 
it is difficult to tell by how much. Strange quark masses as low as ra s (2 GeV) = 80 MeV 



have been reported in the literature |51|], although the errors on unquenched calculations 
are still large. Lacking more precise information on unquenched lattice calculations we 
take as the average lattice value 

ra s (2 GeV) = (110 ± 20) MeV. (2.44) 



which is very close to the one given by Gupta [K 

A large number of determinations of the strange quark mass from QCD sum rules exist 
in the literature. Historically, QCD sum rule results for ra s are given at a scale 1 GeV. 
Taking an average over recent results [p3|-[58 we find ra s (l GeV) = (170±30) MeV. This 



translates to ra s (2GeV) = (124 ± 22) MeV, somewhat higher than the lattice result but 
compatible within the errors. QCD sum rules also allow to derive lower bounds on the 
strange quark mass. It was found that generally ra s (2GeV) > 100 MeV If these 

bounds hold, they would rule out the very low strange mass values found in unquenched 
lattice QCD simulations. 

Finally, one should also mention the very recent determination of the strange mass 
from the hadronic r-spectral function |5"2"|, O which proceeds similarly to the determina- 



tion of a s from r-decays. Normalized at the r mass, the ALEPH collaboration obtains 
m s (m T ) = (176lg?) MeV which translates to ra s (2GeV) = (1701^) MeV. We observe 
that the central value is much larger than the corresponding results from lattice and sum 
rules although the error is still large. In the future, however, improved experimental 
statistics and a better understanding of perturbative QCD corrections should make the 
determination of ra s from the r-spectral function competitive to the other methods. 
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Table 3: The dictionary between the values of m s (fi) in units of MeV. = 340 MeV 
and m c = 1.3 GeV have been used. 



m s (m c ) 


105 


130 


155 


180 


m s (2 GeV) 


90 


111 


132 


154 


m s (l GeV) 


123 


152 


181 


211 



We conclude that the error on m s is still rather large. In our numerical analysis of 
e'/e, where m s is evaluated at the scale m c , we will set 



130 ±25) MeV, 



roughly corresponding to m s (2 GeV) given in (|2.44|) . 



(2.45) 



2.7 Review of B K , B$ ,2) and Bf /2) 
2.7.1 B K 

The renormalization group invariant parameter Bk is defined through 



B K = B K {fj,) la^bi) 



-2/9 



1 + ^r Js 



(2.46) 



(K°\(sd)v-A(sd)v-A\K°) = -B K (fj,)F 2 K m 



K 



[2A7) 



where J3 = 1.895 and J3 = 0.562 in the NDR and HV scheme respectively. 

There is a long history of evaluating Bk in various non-perturbative approaches. The 
status of quenched lattice calculations |64], ^5[ |66| as of 1998 has been reviewed by Gupta 
|[52|| . The most accurate result for B^(2 GeV) using lattice methods has been obtained 
by the JLQCD collaboration 0: B K (2 GeV) = 0.628 ±0.042. A similar result has been 
published by Gupta, Kilcup and Sharpe f55f last year. The APE collaboration []66[ found 
Bk(2 GeV) = 0.66 ± 0.11 which is consistent with 0, [65|| . The final lattice value given 
by Gupta was then 

(Bk) Lattice 

0.86 ±0.06 ±0.06 (2.48) 
where the second error is attributed to quenching. The corresponding result from the 



APE collaboration [66] was Bk = 0.93 ± 0.16. The most recent global analysis of lattice 
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data including also the UKQCD results gives p7 



B K = 0.89 ±0.13 (2.49) 



in good agreement with ( |2.48[ ). 



In the 1/N approach of [[HJ one finds Bk = 0.70 ± 0.10 fl68| , |69| . The most recent 
analysis in this approach with a modified matching procedure and inclusion of higher 



order terms in momenta gives a bigger range 0.4 < Bk < 0.7 [47[ which results from a 
stronger dependence on the matching scale between short and long distance contributions 
than found in previous calculations. It is hoped that inclusion of higher resonances in the 
effective low energy theory will make the dependence weaker. 

QCD sum rules give results around Bk = 0.5—0.6 with errors in the range 0.2—0.3 |70| . 
Still lower values are found using the QCD Hadronic Duality approach (Bk = 0.39 ±0.10) 



7TJ , the SU(3) symmetry and PCAC (Bk = 1/3) J72] or chiral perturbation theory at 
next-to-leading order (B K = 0.42 ± 0.06) J7J. However, as stressed in @ [7|, SU(3) 
breaking effects considerably increase these values. Finally, the analysis in the chiral 



quark model gives a value as high as Bk = 1-1 ± 0.2 [75 



In our numerical analysis presented below we will use 

B K = 0.80 ±0.15 (2.50) 

which is in the ball park of various lattice and large-N estimates. We will, however, discuss 
what happens if values outside this range are used. 

2.7.2 General Comments on B£ /2) and Bf /2) 

As the different methods for the evaluation of these parameters use different values of /t, 
it is useful to say something about their /i- dependence. As seen in ( |2.22| ) and ( |2.30| ) the 
//-dependences of (Qe(A t ))o an d (Qs(a0)2 are governed by the known //-dependence of m s 
and nid and could also in principle be present in B^ 2 \/j,) and B^ 2 \ji). 

Now, as can be demonstrated in the large-N limit, the //-dependence of l/(m s (/i) ± 
md(/t)) 2 in (Qq)o and (Qs)2 is exactly cancelled in the decay amplitude by the diagonal 
evolution (no operator mixing) of the Wilson coefficients ye(fj) and ys{fj) taken in the 
large-N limit. An explicit demonstration of this feature is given in ||. In the large-N 
limit one also finds 

Bi 1/2) = Bf /2) = 1, (Large - N Limit). (2.51) 
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The //-dependence of B^ 2 ^ and B^ 2 ^ for N = 3 and in the presence of mixing with 
other operators has been investigated in [p5fl . This analysis shows that B^ 2 ^ and B^ 2 ^ 
depend only very weakly on /i, when fi > 1 GeV. In such a numerical renormalization 
study the factors B^ 2 ^ and B^ 2 ^ have been set to unity at /i — m c . Subsequently the 
evolution of the matrix elements in the range 1 GeV < /i < 2 GeV has been calculated 

(1/2) (3/2) 

showing that for the NDR scheme Bq and B% were [i independent within an accuracy 
of 2 %. The \i dependence in the HV scheme has been found to be stronger but still below 
6%. Similar weak /i-dependences have been found for B^ 2 ^ and B^ 2 \ 

These findings simplify the comparison of results for B^q 2 ^ and Bj^ obtained by 
different methods. 

2.7.3 B { 6 1/2) and Bf /2) from the Lattice 

The lattice calculations of B^ g and B\ % have been reviewed by Gupta []53J and the 



APE collaboration | 6q| . They are all given at /i = 2 GeV and in the NDR scheme. The 
most reliable results are found for B-f% . The "modern" quenched estimates for these 
parameters are collected in table |] ||52|| . The errors given there are purely statistical. 
The first three calculations use perturbative matching between lattice and continuum, 
the last one uses non-perturbative matching. All three groups agree within perturbative 
matching that Bj 3 ( 2 ^ are suppressed below unity: B^ 2 ^ « 0.6 and B^ 2 ^ 0.8. The 
non-perturbative matching seems to increase these results by about 20%. It is important 
to see whether this feature will be confirmed by other groups. 

Concerning the lattice results for B^( 2 ^ the situation is worse. The old results read 



B^l 2 \2 GeV) = 1.0 ± 0.2 \ 7% [78 1 . More accurate estimates for B^ 2 ^ have been given in 



|7|: B ( i /2) {2 GeV) = 0.67±0.04±0.05 (quenched) and^ 1/2) (2 GeV) = 0.76±0.03±0.05 
(/ = 2). However, as stressed by Gupta |52| , the systematic errors in this analysis are 
not really under control. A recent work of Pekurovsky and Kilcup fj6 |, in which B^ 2 ^ 



is even found to be negative, unfortunately supports this criticism. We have to conclude 
that there are no solid predictions for B^q 2 ^ from the lattice at present. 

2.7.4 B { b 1/2) and Bf /2) from the 1/N Approach 

The 1/N approach to weak hadronic matrix elements was introduced in [0. In this 
approach the 1/N expansion becomes a loop expansion in an effective meson theory. In 
the strict large-N limit only the tree level matrix elements of Qq and Q$ contribute and 
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Table 4: Lattice results for B\ % (2 GeV) obtained by various groups. 



Fermion type 


£(3/2) 


R (3/2) 


Matching 


Staggered 
Wilson [76] 
Clover [ 66 
Clover [IBB] 


0.62(3)(6) 
0.58(2)(7) 
0.58(2) 
0.72(5) 


0.77(4)(4) 
0.81(3)(3) 
0.83(2) 
1.03(3) 


1-loop 
1-loop 
1-loop 
Non-pert. 



one finds ( |2.51| ) while B^ 2 ^ = B^ 2 ^ = 0. The latter fact is not disturbing, however, as 
the operators Q§ and Qi having small Wilson coefficients are unimportant for e'/e. 

(1/2) (3/2) 

In view of the fact that for Bq = B% = 1 and the known value of m t there is 
a strong cancellation between gluon and electroweak penguin contributions to e'/e, it is 
important to investigate whether the 1/N corrections significantly affect this cancellation. 



This has been investigated in [48], where a calculation of (Qe)o and (Qs}2 m the twofold 
expansion in powers of external momenta p, and in 1/N has been presented. The final 
results for (Qe)o and (Qs)2 in ||8] include the orders p 2 and p°/N. For (Q$)2 also the 
term p° contributes. Of particular interest are the O(p /N) contributions resulting from 
non-factorizable chiral loops which are important for the matching between long- and 
short-distance contributions. The cut-off scale A c in these non-factorizable diagrams is 
identified with the QCD renormalization scale fi which enters the Wilson coefficients. 



Table 5: Results for Bq and B$ obtained in the 1/N approach. 





A c = 0.6 GeV 


A c = 0.7 GeV 


A c = 0.8 GeV 


A c = 0.9 GeV 


_g(l/2) 


1.10 


0.96 


0.84 


0.72 




(1.30) 


(1.19) 


(1.09) 


(0.99) 


R (3/2) 


0.64 


0.56 


0.49 


0.42 




(0.71) 


(0.65) 


(0.59) 


(0.53) 



In table |5], taken from |48|, ^7| , we show the values of B^ 2 ^ and B^ 2 ^ as functions of 
the cut-off scale A c . The results depend on whether F n or Fx is used in the calculation, 
the difference being of higher order. The results using Fx are shown in parentheses. The 
decrease of both B-factors with A c = /i is qualitatively consistent with their /i-dependence 



20 



found for \x > 1 in ]2"5| , but it is much stronger. Clearly one could also expect a stronger 



/i-dependence in the analysis of [25[] for /x < 1 GeV, but in view of large perturbative 
corrections for such small scales a meaningful test of the dependence in table [| cannot 
be made. We note that for A c = 0.7 GeV the value of Bq 1 ^ is close to unity as in the 

(3 /2) 

large-N limit. However, B$ is considerably suppressed. An interesting feature of these 
results is the near A c independence of the ratio B^ 2 ^ / B^ 2 \ Consequently the results in 



17 1 can be summarized by 



B 



(1/2) 

- ~ 1.72 (1.84), 0.72 (0.99) < B% ,2) < 1.10 (1.30) . (2.52) 



£(3/2) 



It is difficult to decide which value should be used in the phenomenology of e'/e. On 
the one hand, for A c > 0.7 GeV neglected contributions from vector mesons in the 
loops should be included. On the other hand for A c = fi = 0.6 GeV the short distance 
calculations are questionable. Probably the best thing to do at present is to vary A c = // 
in the full range shown in table [|. This has been done in a recent analysis |50| in which 
e'/e has been found to be a decreasing function of A c . 

Finally, we would like to mention that the first non-trivial 1/N corrections to the 



matrix elements of Qj have been calculated in [81 1 using the methods developed in [82 



(3/2) 

In particular it has been found that B) is a rather strongly increasing function of /x 
with negative values for /i < m c , B\ (m c ) = and positive values for \i > m c . This 
strong /x-dependence of Bj 3 ^ 2 ^ is rather surprising as the numerical renormalization group 
analysis in |25| has shown a rather weak dependence of this parameter. We suspect 
that the inclusion of the full mixing between Q7 and other operators in the analysis of 
| 3T| would weaken the /i-dependence of B\ considerably. While this issue requires an 



additional investigation, the value of B^ 2 ^ has fortunately only a minor impact on e'/e. 
Setting B^ 2 \m c ) = instead of Bf' z '(m c ) = B%\m c ) used here would change our 
results for e'/e only by a few percent. 

2.7.5 B^ 1/2) and Bf ,2) from the Chiral Quark Model 

Effective Quark Models of QCD can be derived in the framework of the extended Nambu- 



Jona-Lasinio model of chiral symmetry breaking ||83|| . For kaon decays and in particular 
for e'/e, an extensive analysis of this model including chiral loops, gluon and 0(p A ) cor- 



rections has been performed over the last years by the Trieste group |84], j3q| . The crucial 



parameters in this approach are a mass parameter M and the condensates (qq) and 



21 



(a s GG). They can be constrained by imposing the AI = 1/2 rule. 



Since there exists a nice review |38j by the Trieste group, we will only quote here their 
estimates of the relevant Bi parameters. They are given in the HV scheme as follows 



B, 



(1/2) 



1.6 ±0.3, 



B 



(3/2) 



0.92 ±0.02 



(Chiral QM). 



Translating these values into the NDR scheme by means of (|2.43|) one finds 



B 



(1/2) 



1.33 ±0.25, 



B. 



(3/2) 



0.77 ±0.02 



(NDR). 



(2.53) 



(2.54) 



We observe a substantial enhancement of B^ 2 " 1 in the chiral quark model, not found in 



other calculations, and a moderate suppression of B 



(3/2) 



The errors given above arise 



from the variation of m B . We will return to this point in subsection [3]5 



It should be remarked that the definitions of the Bi parameters used in |$8| agree 
with our definitions only if in the vacuum insertion formulae in |38[ the (qq) condensate 
is given in terms of m s as follows: 

1 2 



(wr = -f- 



m 



K 



m s ± ma 



(2.55) 



This means that in the usual PCAC relation one has to set Fx = F n . 

It is interesting to observe that in this method B$ / B 8 = 1.74 ± 0.33 in the ball 



park of (2.52). It will be of interest to see whether future lattice calculations will confirm 



this correlation between B, 



(1/2) 



and B. 



(3/2) 



2.7.6 B^ /2) and the AI = 1/2 Rule 



In one of the first estimates of e'/e, Gilman and Wise 1 used the suggestion of Vain- 
shtein, Zakharov and Shifman ]8C| that the amplitude ReAo is dominated by the QCD- 
penguin operator Qq. Estimating (Qe)o in this manner they predicted a large value of 
e'/e. Since then it has been understood |j45| , |46| , |4"?|j that as long as the scale [i is not much 
lower than 1 GeV the amplitude ReAo is dominated by the operators Q\ and Q2, rather 
than by Qq. Indeed, at least in the HV scheme the operator Qq does not contribute to 
Re^o for fi = m c at all, as its coefficient z 6 (m c ) relevant for this amplitude vanishes. Also 
in the NDR scheme z 6 (m c ) is negligible. 

For decreasing /i the coefficient z 6 (/i) increases and the Qq contribution to ReA is 
larger. However, if the analyses in |45|, ^7| are taken into account, the operators Q\ 
and Q2 are responsible for at least 90% of Rev4 if the scale fi = 1 GeV is considered. 
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Therefore in our opinion there is no strict relation between the large value of e' / e and the 
AI = 1/2 rule as sometimes stated in the literature. Moreover, if the 90% contribution 
of the operators Q\ and Q2 to Re^o is taken into account and zq{1 GeV) is calculated 
in the NDR scheme, B^ 2 ^ cannot exceed 1.5 if m s (l GeV) = 150 MeV. Consequently 
we do not think that values of B^ 2 ^ in the NDR scheme as high as 4.0 suggested in 



f40fl are plausible. Unfortunately, due to the very strong fi and renormalization scheme 
dependences of Zq(/i), general definite conclusions about B^ 2 ^ cannot be reached in this 
manner at present. Similarly, we cannot exclude the possibility that B^ 2 ^ is substantially 
higher than unity if it turned out that the present methods overestimate the role of Qi 
and Q2 in Re/Lo- 



2.7.7 Summary 

(3/2') 

We have seen that most non-perturbative approaches discussed above found B\ below 
unity. The suppression of B^ 2 ^ below unity is rather modest (at most 20%) in the lattice 
approaches and in the chiral quark model. In the 1/N approach Bg is rather strongly 
suppressed and can be as low as 0.5. 

Concerning B^ 2 ^ the situation is worse. As we stated above there is no solid prediction 
for this parameter in the lattice approach. On the other hand while the average value of 
.Bg 1 ^ in the 1/N approach is close to 1.0, the chiral quark model gives at fi — 0.8 GeV 
and in the NDR scheme the value for B^ 2 ^ as high as 1.33 ± 0.25. Interestingly both 
approaches give the ratio B^ 2 " 1 / B^ 2 ^ in the ball park of 1.7. 

Guided by the results presented above and biased to some extent by the results from 
the large-N approach and lattice calculations, we will use in our numerical analysis below 
Bi 1/2) and Bf /2) in the ranges: 

B^ /2) = 1.0 ± 0.3, Bf /2) = 0.8 ± 0.2 (2.56) 

(1/2) (3/2) (1/2) 

keeping always B§ > B$ . In our 1996 analysis [0 we have used B y 6 1 ' = 1.0 ± 0.2 
and Bf ,2) = 1.0 ± 0.2 without the constraint B^ /2) > Bf /2) . The decrease of B^ 3/2) 
below unity is motivated by the recent results discussed above. The increase in the range 
of B^ 2 " 1 is supposed to take effectively into account the uncertainty in Q^+rf which we 
estimate to be at most ±30% i.e Q v + V > = 0.25 ± 0.08. We will return to this point in 
Section 3. 
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3 Numerical Results in the Standard Model 



3.1 Input Parameters 

In order to make predictions for e'/e we need the value of ImAt. This can be obtained 
from the standard analysis of the unitarity triangle which uses the data for \V c b\, \V u b\, e, 
AMd and AM S , where the last two measure the size of B® s — B® s mixings. Since this 
analysis is very well known we do not list the relevant formulae here. They can be found 
for instance in || |88|| . 

The input parameters needed to perform the standard analysis of the unitarity triangle 
are given in table |6], where m t refers to the running current top quark mass defined at 
= mf ole . It corresponds to m^ ole = 174.3 ±5.1 GeV measured by CDF and DO [57]. 



We also recall that the lower bound on AM S together with AM^ puts the following 
constraint on the ratio |V^|/|Vt s |: 



\V U \ AMd t _ F Bs ^W s 

\Vu\ ^ m B J AMf^ ~ F Bd JB~B~ d ' ^ 



The range for in table || corresponds roughly to a s (M z ) = 0.119 ± 0.003. 

3.2 Monte Carlo and Scanning Estimates of e'/e 

In what follows we will present two types of numerical analyses of ImA t and e'/e: 

• Method 1: The experimentally measured numbers are used with Gaussian errors 
and for the theoretical input parameters we take a flat distribution in the ranges 
given in table ^. 

• Method 2: Both the experimentally measured numbers and the theoretical input 
parameters are scanned independently within the ranges given in table || 

Using the first method we find the probability density distributions for ImAt and e' / e 
in figs. |T] and ^ respectively. From the distributions in figs. |T] and we deduce the following 
results: 

IrnAj = (1.33 ±0.14) ■ 10" 4 (3.58) 
e'/e = (7.7 ±| g) • 10" 4 (NDR) (3.59) 
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Table 6: Collection of input parameters. We impose ^ 1/2) > B, 



,(3/2) 



Quantity 



Central 



Error 



Reference 



\Vcb\ 

\V ub \ 
B K 
\fB~ d F Bd 
m t 
AM, 
AM S 

MS 

m s (m c ) 



B, 
B 



(1/2) 
6 

(3/2) 



0.040 
3.56- 10- 3 

0.80 
200 MeV 
165 GeV 
0.471 ps" 1 
> 12.4 ps" 1 
1.14 
340 MeV 
130 MeV 
1.0 
0.8 



±0.002 
±0.56 • 10~ 3 
±0.15 
±40 MeV 
±5 GeV 
±0.016 ps- 1 
95%C.L. 

±0.08 
±50 MeV 
±25 MeV 
±0.3 
±0.2 



See Text 

n 



III 



m 



EI 



m 
ii 

See Text 
See Text 
See Text 



Since the probability density in fig. [I] is rather symmetric we give only the mean and 
the standard deviation for lm\ t . On the other hand, the resulting probability density dis- 
tribution for e'/ e is very asymmetric with a very long tail towards large values. Therefore 
we decided to quote the median and the 68% (95%) confidence level intervals. This means 
that 68% of our data can be found inside the corresponding error interval and that 50% 
of our data has smaller e'/ e than our median. 

We observe that negative values of e'/e can be excluded at 95% C.L. For completeness 
we quote the mean and the standard deviation for e'/e: 

e'/e = 9.1 ±6.2 (NDR) (3.60) 

Using the second method and the parameters in table || we find : 

1.04- 10~ 4 < ImAt < 1.63- 10 -4 (3.61) 

1.05 • 10~ 4 < e'/e < 28.8 • 10" 4 (NDR). (3.62) 
The above results for e'/e apply to the NDR scheme, e'/e is generally lower in the HV 

fl/2) ('3/2') 

scheme if the same values for B$ and B$ are used in both schemes. As discussed in 
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2.5 



o 
T3 



2 - 



X 1.5 



o 



1 - 



0.5 - 



ImA t (e) 
lmX t (e') 



1.33 ±0.14) • 10~ 4 
1.38 ±0.14) ■ 10~ 4 



0.4 



0.6 



0.8 




1 1.2 
ImA t [10~ 4 



Figure 1: Probability density distributions for ImA 4 without (solid line) and with (dashed 
line) the e'/e-constraint. 

(1/2) (3/2) 

subsection |2.5| , such treatment of B$ and B% is the proper way of estimating scheme 
dependences at present. 

Using the two error analyses we find respectively: 



e'/e = (5.2 ± 4 : °) ■ 10" 4 (HV) 



and 



0.26 • 10~ 4 < e'/e < 22.0 • 10~ 4 
Moreover, the mean and the standard deviation read 

e'/e = 6.3 ±4.8 (HV). 



(HV). 



(3.63) 



(3.64) 



(3.65) 



The corresponding probability density distribution for e'/e is compared to the one ob- 
tained in the NDR scheme in fig. |^. Assuming, on the other hand, that the values in 
( [2.56j ) correspond to the NDR scheme and using the relation ( |2.43| ), we find for the HV 
scheme the range 0.58 • 10~ 4 < e'/e < 26.9 • 10 -4 which is much closer to the NDR result 
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e'/e [10~ 4 ] 

Figure 2: Probability density distributions for e'/e in NDR and HV schemes. 

in ( |3.62| ). This exercise shows that it is very desirable to have the scheme dependence 
under control. 

We observe that the most probable values for e'/e in the NDR scheme are in the ball 
park of 1 ■ 10 -3 . They are lower by roughly 30% in the HV scheme if the same values for 
(Bq , -Bg ) are used. On the other hand the ranges in (|3.62| ) and ( p.64j ) show that for 
particular choices of the input parameters, values for e'/e as high as (2 — 3) • 10~ 3 cannot 
be excluded at present. Let us study this in more detail. 

3.3 Anatomy of e'/e 
3.3.1 Global Analysis 

In table [7|we show the values of e'/e in units of 10~ 4 for specific values of B^ 2 \ B^ 2 ^ and 
m s (m c ) as calculated in the NDR scheme. The corresponding values in the HV scheme 
are lower as discussed above. The fourth column shows the results for central values 
of all remaining parameters. The comparison of the the fourth and the fifth column 
demonstrates how e'/e is increased when is raised from 340 MeV to 390 MeV. As 

' MS 
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stated in ( |1.7| ) e'/e is roughly proportional to A^. Finally, in the last column maximal 
values of e'/e are given. To this end we have scanned all parameters relevant for the 
analysis of \m.\ t within one standard deviation and have chosen the highest value of 
^Ms = ^® MeV. Comparison of the last two columns demonstrates the impact of the 
increase of lm\ t from its central to its maximal value and of the variation of mt. 

Table ^ gives a good insight in the dependence of e' / e on various parameters which is 
roughly described by (|1.7|) . We observe the following hierarchies: 

• The largest uncertainties reside in m a , b£ I2) and^ 3/2) . e'/e increases universally 
by roughly a factor of 2.3 when m s {m c ) is changed from 155 MeV to 105 MeV. The 
increase of B^ 2 ^ from 1.0 to 1.3 increases e'/e by (55 ±10)%, depending on m s and 
_g( 3 / 2 )_ rj^g correS p 0nc ii n g changes due to B^ 2 ^ are approximately (40 ± 15)%. 

• The combined uncertainty due to ImA t and m t , present both in ImA t and F E t, is 
approximately ±25%. The uncertainty due to m t alone is only ±5%. 

• The uncertainty due to is approximately ±16%. 

The large sensitivity of e'/e to m s has been known since the analyses in the eighties. 
In the context of the KTeV result this issue has been analyzed in ||40||. It has been found 



(1/2) (3/2) 

that provided 2Bq — B% < 2 the consistency of the Standard Model with the KTeV 
result requires the 2a bound m s (2GeV) < 110 MeV. Our analysis is compatible with 
these findings. 

It is of interest to investigate the impact of the relation ( |2.52|) on our results. Scanning 



all parameters in the ranges given in table [] and imposing B^ 2 ^ = 1.7 • B^ 2 ^ we find 

3.7 -10" 4 <e'/e< 26.2 • 10~ 4 (3.66) 



which is somewhat reduced with respect to (|3.62| ) 



Finally we would like to comment on formula ( |1.6| ) in which ReA 4 appears instead of 
ImAt. Since F e decreases with decreasing ReAf one can come closer to the experimental 
data for e' / e by choosing ReAf sufficiently small. In the Wolfenstein parametrization ReA< 
is proportional to 1 — g and a small ReA t corresponds to a sufficiently large positive value 
of the parameter g. Yet it is known from analyses of the unitarity triangle that g is 
bounded from above by the ratio |V^/Vd,| and even stronger by the value of e. If these 
constraints are taken into account the analysis using (|L6| ) reduces to the one presented 
above. 
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Table 7: Values of e'/e in units of 10 4 for specific values of B$ , B% ' and m s (m c ) 
and other parameters as explained in the text. 



£(1/2) 


£(3/2) 


m s (m c )\MeV} 


Central 


AgL = 390 MeV 

MS 


Maximal 






105 


20.2 


23.3 


28.8 


1.3 


0.6 


130 


12.8 


14.8 


18.3 






155 


8.5 


9.9 


12.3 






105 


18.1 


20.8 


26.0 


1.3 


0.8 


130 


11.3 


13.1 


16.4 






155 


7.5 


8.7 


10.9 






105 


15.9 


18.3 


23.2 


1.3 


1.0 


130 


9.9 


11.5 


14.5 






155 


6.5 


7.6 


9.6 






105 


13.7 


15.8 


19.7 


1.0 


0.6 


130 


8.4 


9.8 


12.2 






155 


5.4 


6.4 


7.9 






105 


11.5 


13.3 


16.9 


1.0 


0.8 


130 


7.0 


8.1 


10.4 






155 


4.4 


5.2 


6.6 






105 


9.4 


10.9 


14.1 


1.0 


1.0 


130 


5.5 


6.5 


8.5 






155 


3.3 


4.0 


5.2 



3.3.2 Parametric vs. Hadronic Uncertainties 

One should distinguish between parametric and hadronic uncertainties. Parametric un- 
certainties are related to m t , \V u b\, \V c b\ and A^. One should in principle include m s 
in this list. However, in order to extract m s from the kaon mass one encounters large 
non-perturbative uncertainties. Clearly such uncertainties are also present in the deter- 
mination of \V c b\ and in particular in the determination of IV^I, but they are substantially 

(1/2) (3/2) 

smaller. Hence the hadronic uncertainties discussed below are related to B K , B^*\ Bg'* J 
and m s . 

In table [8] we show ranges for e'/e related to various uncertainties. The parametric 
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uncertainties have been obtained for central values of Bk and m s and two choices of 
(-Bg 1//2 \ -Bg 3/ ' 2 ' ) ). The hadronic uncertainties due to Bk, B^ 2 ^ and B^ 2 " 1 have been found 
by setting all the remaining parameters at their central values. The uncertainty due 
to m s has been shown for two choices of (B$ ,B$ ) and all other parameters set at 
their central values. The last row in table |8] shows the total hadronic uncertainty. It is 
evident from this table that hadronic uncertainties dominate, although the reduction of 
parametric uncertainties is very desirable. 

Table 8: Uncertainties in e'/e in units of 10 -4 as explained in the text. 



Uncertainties 


d(1/2) 


o(3/2) 

n s 


(— / ^)min 


(— / ^jmax 


Parametric 


1.0 


0.8 


5.0 


9.5 


Parametric 


1.3 


0.8 


8.4 


15.1 


Hadronic 






3.0 


13.6 


Hadronic (m s ) 


1.0 


0.8 


4.5 


11.3 


Hadronic (m s ) 


1.3 


0.8 


7.6 


17.9 


Hadronic (full ) 






1.7 


21.3 



3.4 B? /2) -Bi m Plot 

In fig. |3| we show the minimal value of B^ 2 ^ for two choices of m s (m c ) and as a 
function of B^ 2 ^ for which the theoretical value of e'/e is higher than 2.0- 10~ 3 . To obtain 
this plot we have varied all other parameters in the ranges given in table |6|. We show also 
the line corresponding to the relation (|2.52| ). We observe that as long as B^ 2 ^ > 0.6, the 



parameter B^ 2 ^ is required to be larger than unity. This plot should be useful when our 
knowledge of Bg 2 \ B^ 2 \ m s and A^ improves. 

3.5 Approximate Scaling Laws for e'/e 
3.5.1 Preliminaries 

Table contains a lot of information on e'/e. This information can be further extended 
by noting that e'/e depends to a very good approximation on certain combinations of the 
input parameters. This is seen in ( |1.7|) and (12.41 ). Here we want to provide scaling laws 
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(3/2) 



Figure 3: Minimal value of B^ 2 ^ consistent with e'/e > 2.0 • 10" 



1.1 



based on these formulae which allow to obtain from table values for e'/e for different 
sets of input parameters. 



(1/2) R (3/2) 



3.5.2 By, B. 



and 



As seen in (|2.41|), e'/e depends on these important three parameters only through Rq and 



Rs defined in (|2.42j ). Using this property one can for instance immediately find that the 
values for e'/e in the tenth row of table |7| can also be obtained for the set 



B, 



(1/2) 



1.50, 



B. 



(3/2) 



0.90, 



m s (m c ) = 130 MeV . 



(3.67) 



This set of parameters is similar to the input parameters used by the Trieste group 

At this point we would like to remark that in principle the determination of B^ 2 ^ and 
B^ 2 ^ in a given non-perturbative framework could depend on the value of m s . This turns 
out not to be the case in the large-N approach |TjJ |48|, [47]]. In the lattice approach this 
question has still to be investigated. On the other hand there are results in the literature 
showing a strong m<,-dependence of the Bi parameters. This is the case for B^ 2 ^ in the 
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chiral quark model where B^ 2 ^ scales like m s ||38|. Similarly values for B^ 2 ^ calculated 



in 



STJ show a strong m s -dependence. 
In the present paper we have varied (B^ 2 \ B^ 2 ^) and m s independently which is in 
accordance with large-N calculations. This resulted in the following ranges for i? 6 and i?g 

0.5 < Re < 1.95, 0.4 < R$ < 1.5 (3.68) 

which are correlated through their common dependence on m s . 

(1/2) (3/2) 

If (Bg''',Br ') depend on m s these ranges could change. In this context one should 
remark that in the chiral quark model p8f the highest value of Rq corresponds to the 



minimal value of B^ 2 ^ and consequently the comparison of the results from the chiral 
quark model and the large-N approach has to be made with care. 

Finally, it should be remarked that the decomposition of the relevant hadronic matrix 
elements of penguin operators into a product of Bi factors times ~\-/m? s although useful in 
the 1/N approach will become unnecessary in the lattice approach, once matrix elements 
of dimension three will be calculable with improved accuracy. 

3.5.3 Ag, and ImA t 

For a s (M z ) = 0.119 ± 0.003, the ratio e'/e is within a few percent proportional to A™ 

JV1 o 

On the other hand e'/e is exactly proportional to ImAt at fixed m t . However, if m t is 
varied the correlation in m t between ImAt extracted from e and F e > has to be taken into 
account. Consequently the simple rescaling of e'/e with the values of ImAf is only true 
within a few percent. 

3.5.4 Sensitivity to Vl^+rf 

The dependence of e'/e on can be studied numerically by using the formula ( |2.4| ) 

or incorporated approximately into the analytic formula fl2.38| ) by simply replacing B^ 2 ^ 
with an effective parameter 

f R (l/2h _ p(l/2) (l -0-9 gg+gQ , s 

(B 6 j cfr - B Q (3.o9J 

A numerical analysis shows that using ( 1 — Q v +r]' ) overestimates the role of Vt v+V i . In our 
numerical analysis we have incorporated the uncertainty in Q v +rj> by increasing the error 
in B { b 1/2) from ±0.2 to ±0.3. 
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The last estimates of Vt^+rf have been done more than ten years ago [IB|-[I8| and it is 
desirable to update these analyses which can be summarized by 



n v+r]/ = 0.25 ± 0.08 . (3.70) 

The uncertainty in e'/e due to fl v+v ' alone is approximately ±12% and is slightly lower 
than the one originating from A^. 

3.5.5 Sensitivity to B K 

As ImA 4 extracted from e increases with decreasing Bk, there is a possibility of increasing 
e'/ e by decreasing Bk below the range considered in table ||. It should be remarked that 
e'/e is not simply proportional to 1/Bk as the extraction of lmX t from e involves also 
KeXt (see (|1.2| )). For the phase 5 in the first quadrant as favoured by the analyses of the 
unitarity triangle fl8~8"| , the dependence of e'/e on B K is weaker than l/B K [p5| . 



Now, the highest value of lm\ t consistent with the unitarity of the CKM matrix is 
1.73 • 10~ 4 . It is obtained from e for B^ = 0.52. This increase of ImAj beyond the range in 
( |3.61| ) would increase the maximal values in table |7| by approximately 6%. On the other 



hand it should be emphasized that for Bk = 0.8 — 0.9, as indicated by lattice calculations, 
e'/e is generally smaller than found in our paper unless B^ 2 ^ is substantially increased. 



This is what happens in the chiral quark model [ 38) where on the one hand B^ = 1.1 ±0.2 



and on the other hand b9~' 2 ^ = 1.6 ± 0.3. 



3.6 Impact on lm\ t and the Unitarity Triangle 

As we stressed at the beginning of this paper the main new parameter to be fitted by 
means of e'/e is ImA*. Our analysis indicates that the Standard Model estimates of e'/e 
are generally below the data. If the parameters m s , B^ 2 \ B^ 2 \ A^ and Q v +rj' are such 



that lm\ t consistent with e (see ( |3.61| )) cannot accomodate the experimental value of e'/e, 
one has to conclude that new contributions from new physics are required. On the other 
hand if the data on e'/e can be reproduced within the Standard Model, then generally 



a lower bound on lmX t excluding a large fraction of the range (|3.61| ) can be obtained. 



Unfortunately, the strong dependence of the lower bound on the parameters involved 
precludes any firm conclusions. Similar comments apply to the possible impact of e'/ e on 
the analysis of the unitarity triangle: the presently allowed area in the (g, fj) plane [Q 
can be totally removed or an improved lower limit on fj from e' / e will decrease the allowed 
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region considerably. As an illustration we show in table |5] the lower bound on Im\ t from 
e'/e as a function of B^ 2 ^ for m s (m c ) = 105 MeV, B^ 2 ^ 



1.3 and A§ 



390 MeV. 



To this end we have used the formula ( |2.2j ) with e'/e > 2.0 • 10~ 3 . Comparing with ( |3.61|) 
we indeed observe that the lower bound on ImAj has been improved. 

The impact of e'/e-data as given in ( |1.9[ ) on lm\ t can also be investigated by the 
method 1 which was used to obtain (|3~58p and d3759|) . We find 



ImAt = (1.38 ±0.14) ■ 10" 



(3.71) 



The corresponding distribution is compared with the one without e'/ ^-constraint in fig. [l|. 
We observe a very modest but visible shift towards higher values for ImA 4 . 



Minimal values of lra\ t , Br(Ki 

(1/2) _ n o i a (4) 



Table 9: 

m s (m c ) = 105 MeV, B 6 
assuming e'/e > 2.0 ■ 10~ 3 



-k°vv) and Br(Ki 



n°e + e ) dir 



1.3 and A^ 



390 MeV and specific values of B$ 



for 

(3/2) 



Bf /2) 


(ImA t ) min 


Br{K L -> 7r°z/z/) min 


Br{K L -> vr°e + e-)^ n 


0.6 
0.8 
1.0 


1.14- KT 4 
1.27- 10" 4 
1.42 • KT 4 


1.8 • 10~ n 
2.2 • 10" 11 
2.7- 10~ n 


3.0 • 10~ 12 
3.7- 10" 12 
4.7- 10~ 12 



3.7 Impact on Kl —> irvv and Ki —> 7T e + e 

The rare decay Kl — > ts^vv is the cleanest decay in the field of K-decays. It proceeds 
almost entirely through direct CP violation |94| and after the inclusion of NLO QCD 
corrections |)5] the theoretical uncertainties in the branching ratio are at the level of 
1 — 2%. Similarly the contribution of direct CP-violation to the decay Kl — * 7c°e + e~ is 
very clean. Using the known formulae for these decays || |)5[ and scanning the parameters 
given in table |] we find: 

1.6 • 10~ n < Br(K L -> ttW) < 3.9 ■ 10~ n (3.72) 

2.8 ■ 10~ 12 < Br{K L -> 7r e + e -) dir < 6.5 ■ 10~ 12 (3.73) 

Since these branching ratios are proportional to (ImA t ) 2 any impact of e'/e on the latter 
CKM factor will also modify these estimates. We illustrate this in table |9| where an 



34 



improved lower bound on ImAt implies improved lower bounds on the branching ratios 
in question. With decreasing B^ 2 ^ and increasing m s these lower bounds continue to 
improve excluding a large fraction of the ranges in ( 3.72|) and ( p.73|) . In obtaining the 
results in table |9| correlations in m t and the CKM parameters between e'/e, e and the 
branching ratios for the decays considered have been taken into account. Unfortunately, 
due to large hadronic uncertainties in e'/e, no strong conclusions can be reached at present. 

In the future the situation will be reversed. As pointed out in p6f the cleanest mea- 
surement of ImAt is offered by Br^K^ 



7T l>l>) 



ImXt = 1.41 • 10" 



4 


"165 GeV" 


1.15 


~Br(K L -> tt°vv) 


1/2 




_ m t {m t ) _ 




3- 10- 11 





(3.74) 



Once ImAt is extracted in this manner it can be used in e'/e thereby somewhat reducing 
the uncertainties in the estimate of this ratio. 



4 Implications for Physics Beyond the Standard Model 
4.1 General Comments 

We have seen that the Standard Model estimates of e'/e are generally below the experi- 
mental results from NA31 and KTeV. In view of the large theoretical uncertainties it is, 
however, impossible at present to conclude that new physics is signaled by the e'/e-data. 
Still, we can make a few general comments on the extensions of the Standard Model with 
respect to e'/e: 

• In models where the phase of the CKM matrix is the only source of CP viola- 
tion, the modifications with respect to the Standard Model come through new loop 
contributions to e and e'/e. If the new contributions to e are positive and the 
contributions to F e i are negative, then ImAj, F e i and consequently e'/e are smaller 
than in the Standard Model putting these models into difficulties. An example of 
this disfavoured situation is the two-Higgs doublet model II in which e'/e has been 



analysed a long time ago [0. We will update this analysis below. 

In the Minimal Supersymmetric Standard Model, the last analysis of e' / e after the 



top quark discovery has been performed in p8[ . Here in addition to charged Higgs 
exchanges in loop diagrams, also charginos contribute. The chargino contribution 
to e has always the effect of decreasing ImA t . However, depending on the choice 
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of the supersymmetric parameters, the chargino contribution to F e i can have either 
sign. Consequently, e'/e in the MSSM can be enhanced with respect to the Standard 
Model expectations for a suitable choice of parameters, low values of chargino (stop) 
masses and high charged Higgs masses. Yet, as stressed in |9*8| , generally F E > is 



further suppressed by chargino contributions and the most conspicuous effect of 
minimal supersymmetry is a depletion of e'/e. 

The situation can be different in more general models in which there are more 
parameters than in the two Higgs doublet model II and in the MSSM, in particular 
new CP violating phases. As an example, in more general supersymmetric models 



e'/e can be made consistent with experimental findings [[42], p9[ . Unfortunately, in 
view of the large number of free parameters such models are not very predictive. 
Similar comments apply to models with anomalous gauge couplings [f|l[] and models 



with additional fermions and gauge bosons ||100|| in which new positive contributions 
to e' / e are in principle possible. A recent discussion of new physics effects in e' j e can 
also be found in ||40| . In the past, there have of course been several other analyses 
of e'/e in the extensions of the Standard Model but a review of these analyses is 
clearly beyond the scope of this paper. 



Finally, models with an enhanced sdZ vertex, considered in ||101||, can give rise to 



large contributions to e'/e as pointed out in ||102 . As analyzed in the latter paper 



in these models there exist interesting connections between e'/e and rare K decays. 



4.2 An Update on e'/e in the Two-Higgs Doublet Model II 



A detailed renormalization group analysis oie'/e in the Two-Higgs Doublet Model II 
has been presented in |JTJ . It has been found that due to additional positive charged Higgs 
contributions to e and corresponding negative contributions to F e i through the increase 
of the importance of Z°-penguin diagrams, the ratio e'/e is suppressed with respect to 
the Standard Model expectations. Since this analysis goes back to 1990 and several input 
parameters, in particular m t , have been modified we would like to update this analysis. 

We recall that the two new parameters relevant for our analysis are the charged Higgs 
mass (Mh) and tan/3, the ratio of the two vacuum expectation values. The expressions 
for the new contributions with charged Higgs exchanges to e are rather complicated and 



will not be repeated here. They can be found in Section 3 of [f)7|j . The QCD corrections to 
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these contributions are given there in the leading logarithmic approximation. As of 1999 
only NLO corrections to box diagram contributions with internal top-quark exchanges 
to B° — B° mixing are known p8 |. Unfortunately the NLO QCD analysis for e in the 
2HDMII is still lacking. For this reason we have used the leading order expressions for the 
Higgs contributions to e |57| except for the box diagram contributions with internal top- 
quark exchanges where we took the NLO QCD factor obtained in the Standard Model. 
While such a treatment is clearly an approximation, it is sufficient for our purposes. 

The analysis of F £ i on the other hand can be done fully at the NLO level. We only 
have to add to the functions X Q (x t ), Y (x t ), Z (x t ) and E (x t ) the contributions from 
charged Higgs exchanges. They are given as follows: 



AX = AY 



tan 2 j3 



}(y-l) 8(x-l) : 



logy 



where 



D H (y) 



AZn = AXn 



1 



4 tan 2 (3 
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y{3y 



tan 2 /3 

y (47y 2 
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+ 



y 



6(y- 
-3y 3 + 6y-A) 



logy 



logy 



(4.75) 

(4.76) 
(4.77) 

(4.78) 



108 (y - 1)" 18 (y - l) 4 

with y = m 2 /M H 2 . 

We observe that all new contributions to F e > are inversely proportional to tan 2 f3. In e 
they are inversely proportional to tan 2 (3 and tan 4 (3. This should be contrasted with the 
case of B — > where there are new contributions with charged Higgs exchanges, which 
do not involve tan/5. Thus e'/e is more sensitive to tan/? than B — > X s 7. This implies 
that in principle a better constraint for tan (3 could be obtained from e' / e than from the 
latter decay. 

It is obvious from this discussion that e' / e in the 2HDMII is lower than in the Standard 
Model for any choice of input parameters. Consequently, for low Mh and tan (3, the ratio 
e'/e is generally well below the experimental data. On the other hand if the Standard 
Model is consistent with the experimental value of e'/e, it is possible to put a lower 
bound on tan/? as a function of Mh- In fig. f| we show the result of such an analysis 
for B^ 2 ^ = 0.8 and selected values of B^ 2 ^ and m s . The remaining parameters have 
been scanned in the ranges given in table |6|. We require e'/e > 2.0 ■ 10~ 3 . We observe 
that for the lowest values of M H « 200 GeV allowed by the B -> X^ decay [p3j -p5|, 
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Figure 4: Lower bound on tan (3 as a function of Mh consistent with e'/e > 2.0 ■ 10 



m s (m c ) = 105 MeV and Bq = 1.3 the lower bound on tan/? is similar to the one 
obtained from B — > X s 7. For higher values of m s (m c ) and lower values of B^ 2 ^ the 
bound on tan /3 becomes stronger than from — > X s 7. 



5 Conclusions and Outlook 

We have presented a new analysis of e'/e in the Standard Model in view of the recent 
KTeV measurement of this ratio, which together with the previous NA31 result firmly 
establishes direct CP violation in nature. Compared with our 1996 analysis ||34)| , the 
present analysis uses improved values of \V u b\, \V c b\, m t , and m s as well as new 
insights in the hadronic parameters B^ 2 \ B^ 2 ^ and Bk- Our findings are as follows: 

• The estimates of e'/e in the Standard Model are typically below the experimental 
data. Our Monte Carlo analysis gives 

(NDR) 




e'/e={ _ ) m J (5-79) 

-2.7) 



4-6^. 10 -4 ( HV ) 
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The difference between these two results indicates the left over renormalization 
scheme dependence. 

• On the other hand a simple scanning of all input parameters gives 

1.05 • 1(T 4 < e'/e < 28.8 • 10" 4 (NDR) (5.80) 

and 

0.26 • 10~ 4 < e'/e < 22.0 • 10~ 4 (HV) (5.81) 

This means that for suitably chosen parameters, e' / e in the Standard Model can be 
made consistent with data. However, this happens only if all relevant parameters 
are simultaneously close to their extreme values. This is clearly seen in table |7| and 
fig. |3|. Moreover, the probability density distributions for e' / e in fig. ^| indicates that 
values of e'/e in the ball park of NA31 and KTeV results are rather improbable. 

• Unfortunately, in view of very large hadronic and substantial parametric uncertain- 
ties, it is impossible to conclude at present whether new physics contributions are 
indeed required to fit the data. Similarly it is difficult to conclude what is precisely 
the impact of the e' / e-data on the CKM matrix. However, there are indications as 
seen in table ^| that the lower limit on ImAt is improved. The same applies to the 
lower limits for the branching ratios for — > tt°uj? and — > 7r°e + e~ decays. 

• Finally, we have pointed out that the e' / e data puts models in which there are new 
positive contributions to e and negative contibutions to e' in serious difficulties. In 
particular we have analyzed e' / e in the 2HDMII demonstrating that with improved 
hadronic matrix elements this model can either be ruled out or a powerful lower 
bound on tan /3 can be obtained from e'/ e. 

The fact that one cannot firmly conclude at present that the data for e'/e requires 
new physics is rather unfortunate. In an analogous situation in the very clean rare decays 
K — ► -nvv a departure of the experimental result from the Standard Model expectations 
by only 30% would give a clear signal for new physics. This will indeed be the case 
if the improved measurements of the Br(K + — > n + uu) from BNL787 collaboration at 
Brookhaven | |106| find this branching ratio above 1.5 ■ 10~ 10 . All efforts should be made to 
measure this branching ratio and the branching ratio for Kl — > 7r°z/j/, which while being 
directly CP violating is almost free of theoretical uncertainties. 
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The future of e' / e in the Standard Model and in its extensions depends on the progress 
in the reduction of parametric and hadronic uncertainties. We have analyzed these un- 
certainties in detail in Section 3 with the results given in table |8|. 

Concerning parametric uncertainties related to \V c b\, m t and A^, we expect 

that they should be reduced considerably in the coming years. This will, however, result 
only in a modest reduction of the total uncertainty in e' / e. In this respect a measurement 
of ImA t in a very clean decay like Kl — ► 7r°z/z/ would be very useful. 

A real progress in estimating e' / e will only be made if the non-perturbative parameters 
B K , Bg 2 \ Sg 3//2 ' ) and £1^+^ as well as the strange quark mass m s will be brought under 
control. In particular the sensitivity of non-perturbative methods to /i and renormaliza- 

(1/2) (3/2) 

tion scheme dependences of B & and B 8 is clearly desirable. We expect that consid- 

A (3/2) 

erable progress on B^ and B$ should be made in the coming years through improved 
lattice calculations. Progress on f^+^z should also be possible in the near future. More- 

(3/2) 

over, as various estimates of Bk, B$ and £l v+v > by means of several non-perturbative 
methods are compatible with each other we do not expect big surprises here. Similar 
comments apply to \V u b\, \V c b\, m t and A™.. 

On the other hand, it appears that it will take longer to obtain acceptable val- 
ues for m s and B^ 2 \ In view of the bounds p9|-||61||, it is difficult to imagine that 
^s(^c) < 105 MeV. Consequently we expect that future improved estimates of m s will 
most probably exclude the lowest values of m s considered in this paper. This would si- 
multaneously exclude the highest values for e'/e obtained by us unless B^ 2 ^ is found to 
be higher than used here. In this respect improved estimates of B^ 2 \ if found substan- 
tially higher than unity, could have considerable impact on our analysis. Finally, it should 
be stressed that future lattice calculations will give the full matrix elements without the 
necessity to use separately (B^ 2 \ B^ 2 ^) and m s . 

In any case e'/e already played a decisive role in establishing direct CP violation in 
nature and its rather large value gives additional strong motivation for searching for this 
phenomenon in cleaner K decays like Ki — > 7c uu and — > 7r°e + e~, in B decays, in D 
decays and elsewhere. 
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